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We re-examine one-pion exchange as a possible binding mechanism in the X(3872) charmonium-like state 
and argue it to be not sufficiently binding for this purpose. We conclude therefore that other short-range 
dynamics are responsible for the X formation. 
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1. INTRODUCTION 

A charmonium-like state X(3872) was observed by 
Belle Collaboration in 2003 [1 in _B-meson decays in 
the mode B + — > K + X — > K + ir + ir~ J/ip, and since 
then this state attracts attention of both experimental- 
ists and theorists due to very specific properties it pos- 
sesses. This state is also seen in the mode 7r + 7r~7r° J/ip 
[2], with approximately the same branching fraction as 
that in the mode tt +, k~ J/i/j (for a recent update see 
[3]), that clearly indicates a large isospin violation in 
the X. There is no general agreement on the quan- 
tum numbers of this state yet: while the analysis of 
|the Tr + ir~J/ip decay mode yields either J PC = 1 ++ or 
J PC = 2 h quantum numbers [3111], the recent anal- 
ysis of the ir + Tr~ir°J/ip mode seems to favour the 2 h 
assignment ^2J, though the 1 ++ option is not excluded. 
Quite opposite claims can be found in the literature (see, 
for example, [5] versus [5]), however there is a general 
consensus that, if 2 ^ quantum numbers are confirmed, 
Very exotic theoretical models for the X will have to be 
invoked to explain all established properties of the lat- 
ter [2l 13 ffl. In addition, A (3872) resides exactly at 
the D°_D* c lij threshold that, for the quantum numbers 
1 ++ , can be explained naturally by a strong attraction 
of a conventional 2 3 Pi cc charmonium to the S'-wave 
DD* threshold QUI EE]. Molecule model [El Q31 M CGO 
is also consistent only with the option 1 ++ . In this 
Letter we stick to the quantum numbers 1 ++ , as to the 
most promising candidate, and revisit the long standing 
problem of the one-pion exchange (OPE) as a possible 
binding mechanism responsible for the formation of the 
X. 



1 'A proper C-parity eigenstate is always meant by this (and 
similar) shorthand notation. 



Pion exchange between charmed mesons was sug- 
gested long ago [T2l [13] as a mechanism able to bind the 
isosinglet DD* mesonic system and to form a deuteron- 
like state near threshold. This model was revisited 
shortly after the A(3872) discovery [HUS], while fur- 
ther implications of the nearby pion threshold are dis- 
cussed in [THl E]. According to a combined analysis 
of the most recent data for the X —> tt + tt~J/i/j and 
X -> D a D a n a modes HB1 Q33 H3] , A(3872) is a bound 
state localised within approximately 1 MeV below the 
neutral DD* threshold. Therefore a relevant question 
arises concerning the mechanisms responsible for the 
formation of this bound state, and the OPE is a possi- 
ble candidate for this role. This problem is addressed 
in [2T1[22]. In particular, only the neutral D°D*° config- 
uration is studied in |21] and the conclusion made is in 
negative, namely that the OPE appears to be unable to 
bind such a system. However, in [22], the charged DD* 
channel is also taken into account and the opposite con- 
clusion is made that the resulting OPE interaction is 
able to produce the X as a shallow bound state. Notice 
that both above-mentioned calculations treat the DD* 
system in a deuteron-like fashion: OPE enters in the 
form of a static potential. In the meantime, the D*° 
mass is very close to the D°tt° threshold and thus a rel- 
evant worry [53] is that in the DD* system, bound by 
the OPE, the pion may go on shell. The latter calls for 
the proper inclusion of the three-body DDn unitarity 
cuts. As shown in [23], the cut effects are of paramount 
importance in the D a Dp system, if the Dp width is 
dominated by the S'-wave Dp — » D a n decay. In par- 
ticular, (deeply) bound D\D* states found in [33 [2S] 
in the static approximation disappear completely from 
the spectrum if the full three-body treatment is invoked. 
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In the case of the AT (3872) one deals with the P-wave 
D* — > Dtt vertex, so no disastrous consequences are ex- 
pected from the inclusion of the cut effects (see [27J for 
a detailed discussion of the role played by the three- 
body dynamics in the X(3872)), however the problem 
has to be treated with lots of caution because of the 
divergent D-meson loop integrals. In [2T] [22] divergent 
integrals are made finite with the help of form factors of 
a suitable form, with a cut-off parameter A, and the con- 
clusion whether the OPE interaction is binding enough 
is made as based on the value of A necessary to pro- 
duce a bound state with the given binding energy. In 
particular, it is shown in [21] that a bound state in the 
D°D*° system, with the binding energy around 1 McV, 
exists only for the values of A of order of 6 GeV, that is 
for values much larger than those admitted by interpre- 
tation of the form factors in terms of the quark model. 
On the contrary, it is argued in [22] that, already for as 
small cut-off's as just 1 GeV, a bound state with the 
binding energy 1 MeV appears, if the charged channel 
is taken into account (and provides an extra attraction 
in the system). This result is interpreted then in [22] 
as a proof that the OPE provides enough attraction to 
produce a bound state. 

In our research we also adopt the above mentioned 
regularisation scheme, that is we introduce a form factor 
with the cut-off A, interpret values of A below 1 GeV as 
phenomenologically adequate and disregard larger val- 
ues of A as unphysical. Indeed larger values of the cut- 
off cannot be justified in the framework of quark mod- 
els. Then, in our model we (i) include both neutral and 
charged DD* channels, (ii) go beyond the static approx- 
imation used in [211 122) and include the imaginary part 
of the OPE into consideration, (iii) treat pions as rel- 
ativistic particles, but stick to nonrelativistic dynamics 
for the D and D* mesons, as prompted by the values of 
the loop momentum of order of the cut-off A < 1 GeV. 

The result of this improved calculation which we re- 
port here is that no bound state exists in the system for 
A ~ 1 GeV, and it appears at threshold for as large val- 
ues of the cut-off as A > 2 GeV. Thus the conclusion we 
are led to is that OPE apparently fails to be sufficiently 
binding to produce the A" as a bound state. Therefore 
different, short-range, mechanisms are responsible for 
the X formation. 



2. DD* OPE POTENTIAL 



where is the pion 4-momentum, u* , u, and ir are 
the isospin wave functions of the D*, D 7 and the pion, 
respectively. Then various D* — > Dtt decay widths can 
be evaluated as 



r{D*° 
r{D*+ 
r(D* + 



D°TT°) = 



D+7T ) = 



fl/gpo 
247rmJ ' 

247rm2 c ' 
127rmL ' 



(2) 
(3) 
(4) 



where qoo, q c o, and qo c are the D°-ir°, D + -ir Q , and D°- 
7T + relative momentum, respectively. Notice that, up 
to a kinematical factor, r(D* + — > D°ir + ) is two times 
larger than P(P*° -> D°Tr°) and r(D*+ -> D+tt ). 

Here and in what follows too, m c , m*p, m* c , m^o, 
and m-n-c are the masses of the neutral and charged D 
mesons, D* mesons, and pions, respectively. In our cal- 
culations we use the following values [28] : 

m^o = 134.98 MeV, = 139.57 MeV, 

too = 1864.84 MeV, m c = 1869.62 MeV, 
to»o = 2006.97 MeV, m» c = 2010.27 MeV. 

To have a better contact with previous works, we 
introduce an effective coupling parameter Vq such that 
(see [22] [25]) 



r(D*+ -> D°7r+) = 2V ^, 
mi. 



r{D*° 

and, therefore, 



loo 



A 

Ami 



671% 



(5) 
(6) 

(7) 



where the mass difference between the charged and neu- 
tral states is neglected in the denominators. 

The effective coupling Vo can be estimated from the 
data on the total D*+ width [25], 



r c = (96 ± 22) keV, 

and the branching fractions [25] 

B(D*+ -> D+n ) = (30.7 ± 0.5)%, 
B(D*+ -+ D°tt + ) = (67.7 ± 0. 



(8) 



(9) 
(10) 



We start from the discussion of the P-wave D* 
Dtt vertex which can be parametrised in the form 



(1) 



The value V = 1.3 MeV [22] complies well with the 
data on the D* + pionic decays. 

The total D*° width is not known experimentally, 
however it can be estimated from ([2J, ([3]), and with 
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the help of the branching fraction B(D*° 
(61.9 ±2.9)% [28] to be 

r Q = (65 ± 15) keV. 



D tt ) = 



(11) 



Summarising the discussion above, we can build the 
generic D*Dtt vertex in the form 



g(q) = gq 



A 2 



A 2 



g 



3/2 ' 



(12) 



where q is the relative momentum in the Dir system, 
and a form factor with the cut-off parameter A is intro- 
duced, as was explained above. 

The OPE potential V in the coupled (D*D)-(D*D) 
system can now be built either in the covariant formal- 
ism for the propagator of the intermediate pion or in 
the framework of the Time-Ordered Perturbation The- 
ory. In the latter case it is given by the sum of two 
possible orderings, 



V ^V 1 +V 2 , 



(13) 



which correspond to the DDtt and D*D*tt intermedi- 
ate state, respectively (see Fig. [1} . In what follows we 
consider the C-even isosinglet channel. 

The static limit of the potential (fT3"|) can be obtained 
in the standard manner. If one neglects isospin breaking 
due to the mass difference between charged and neutral 
mesons, then the covariant form of the static po- 
tential in the momentum space is 



VX (<z) 



gn(q)g n >(q) 



(14) 



(27r) 3 q 2 + [m 2 — (to — to*) 2 ] ! 

where the factor 3 accounts for the sum of the neutral- 
and charged-pion exchanges, the spin indices n and n' 
are to be contracted with polarisation vectors of the 
D* mesons, to and m* are the masses of the D and 
D* mesons, respectively (charged and neutral masses 
being indistinguishable in the exact isospin-conserving 
limit). Notice that the potential (|T4"]) is two times 
smaller than the one used in 22 , so that, effectively, the 
coupling parameter Vq used in [35] is two times larger 
than Vq — 1.3 MeV compatible with the data on the 
D* meson pionic decays. Thus we are forced to use the 
value Vq — 2.6 MeV to reproduce the results of [2"2"] . 

In the framework of the Time-Ordered Perturbation 
Theory expression (|T4"]) can be presented as 

Kta"'(9) = - 7 ^9n(q)9n'(q)(vr t (q) + vrHq) 



(2tt) 2 
Vf at (q) 



V 2 stat (q) 



1 



2E 7I (E 7r + to — to*) ' 
1 

2Ei T (E 7! + to, — to) ' 



(15) 
(16) 




Figure 1. Diagrams in the Time-Ordered Perturbation 
Theory responsible for the potentials Vi (left plot) and 
Vi (right plot) from (TO)) . The double and single solid 
line are for the D* and D, respectively, while the dashed 
line is for the pion. The thin vertical line pin-points the 
intermediate state. 



where E^ — y/q 2 + to 2 . 

Since to* rs m + m v , then one has V% <C V\ for the 
momenta |q| of order a few hundred MeV. On the other 
hand, as seen from expression (|15[) . the static approxi- 
mation for the potential V\ is inadequate, and the full 
treatment of the DDtt intermediate state is required. 



3. FORMULATION OF THE PROBLEM 

In [27] the role of the dynamical pions played in 
a near-threshold resonance (at the example of the 
AT (3872) charmonium) was discussed in detail. The key 
idea of the formalism used in [37] was to separate the 
short- and long-range interaction and to study in de- 
tail the role played by the long-range part of the OPE. 
Meanwhile the short-range part of the OPE (together 
with other possible short-range interactions, such as, 
for example, s-channel coupling to charmonium) was 
absorbed into a constant interaction described by a 
counter-term. For a given value of the cut-off A, the 
counter-term was tuned to guarantee the existence in 
the system of a bound state with the given binding en- 
ergy. As a result, the leading dependence of physical ob- 
servables on A was absorbed by the counter-term. This 
approach, standard for effective theories, allowed one to 
use values of the cut-off as small as just a few hundred 
MeV and therefore to deal with nonrelativistic kinemat- 
ics for all particles involved and to neglect the second 
ordering term as well as all contributions coming from 
higher multiparticle intermediate states. On the con- 
trary, the aim of the present research is to investigate 
the OPE (its short-range part) from the point of view of 
the possibility to bind the DD* system alone. To this 
end, as was explained above, we consider the values of 
A < 1 GeV, as prompted by quark models. This implies 
that, while D mesons can still be treated nonrelativis- 
tically, we are forced to resort to relativistic kinematics 
for the pions. Technically this is the only difference be- 
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Figure 2. Kinematics of the DD* scattering due to the 
OPE. Double lines denote D* 's while single lines denote 
D's. 



tween the formulae used in this Letter and those found 
in [37], while their derivation remains the same. We 
therefore skip here all details of the formalism (they can 
be found in [27]) and only pin-point the essential differ- 
ences between the two approaches. In particular, the 
OPE interaction graph depicted in Fig. [2] is described 
by the potential 



1 9n ip' + ot ik p)g n i (p + a ik p') 



(2tt) 3 



D 3 ik(p,p') 



(17) 

where the vertex function g(q) is given in (|12D above, 
while n, n' are i, k are spin and channel indices, respec- 
tively. The three-body propagator reads 



D 3 ( P) p') = 2E„{E„- n-iD), 



where 



fi = m*o + % + E — \J to 2 + p 2 — \/ m' 2 + p' 2 



and 



E n = y/(p + p') 2 + m 2 n . 



(18) 



(19) 



(20) 



The energy E is defined with respect to the D°D* a 
threshold. The coefficients a and a' can be extracted 
from the standard relativistic textbook formula for 
the relative momentum written in terms of the single- 
particle momenta in the form (in notations of Fig. [2] 
see, for example, [501 I5T] ): 



1 



\J e' 2 — p 2 
1 

\J e 2 — p' 2 



pp 



where 



e = \p< 



sj m 2 



E v , 



e' + yje' 2 - p 2 



PP 



e' = yjm' 2 +p> 2 + E„. 



(21) 
(22) 

(23) 



J , , , , 1 , , , , 1 , , , , 1 , , , , 1 p [MeV] 

500 1000 1500 2000 

Figure 3. Dependence of the coefficient a on the mo- 
mentum p (p' = p). 



picted in Fig. [5] read: 
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m , to* 



: TO^O, 



too, m — to c , m^ — m^a^ 



TO C , TO = TOo, to* 

where the channels are defined as: 



10) 



D°D*°, 
D+D*- 



|0> 



D°D*°, 
= D~D*+. 



The potentials V™ are to be supplied by the isospin 
factors: 

(0|r 1 -r2|0)=<c|ri-f 2 |c) = l, (24) 
for the neutral pion exchange, and 



\n-T 2 \c) = (c\n -f 2 |0) = 2, 



(25) 



Nonvanishing components of the OPE potential de- 



fer the charged pion exchange, which we take into ac- 
count explicitly in the Lippmann-Schwinger equations 
(see the system (|2l)|) below). 

The static approximation is obtained from expres- 
sion (|17|) by setting a = a' — 1 and [i = to*o + too — 
to — to'. In Fig. [3J we plot the dependence of the coeffi- 
cient a on the momentum p. It is clearly seen from this 
plot that a decreases fast with the increase of the mo- 
mentum p that, by virtue of (fl~2"|) and (fT7|) . leads to an 
effective suppression of the OPE potential as compared 
to the static approximation with a — a' = 1 for all 
momenta. In this approximation, the potential is local, 
that is, it depends on q = p + p' and, in the isospin- 
conserving limit, the first ordering static potential (|T5|) 
is immediately reproduced. 

Then we can write the system of coupled Lippmann- 
Schwinger equations for the DD* ^-matrix elements 

aoo'(P>P') and a?o"'(P:P') ( see E2 for the details). In 
the C-even channel it reads 
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f a™(p,p',E) =V T'(P,P') 



d 3 ! 



-V T{p,s)a^ n \s,p',E) 



2 ' ^ ) Kc m {P^)a^'{s,p' : E) 



a% (p,p\E) =2V C T (P,P' 
d 3 s 



(26) 



A (s 
d 3 s 



V£r(p,s)a™ n (s,p',E) 
V™(p,s)a2 n '(s,p',E). 



A c (s) 

The inverse propagators Aq and A c take the form 



(27) 



where S = m* c +m c — 77j*o~ m o = 8.08 MeV, the reduced 
masses are defined as 

m*omo m» c m c 



Ho 



t'c 



m*o + tuo ' ' m* c + m c ' 

and the widths To and r c are given in (fTTj) and ©, 
respectively. 

Notice that both two-body propagators Ao(p) and 
A c (p) from ((27|) as well as the three-body propagator 
D^{p : p') from (fT5|) generate contributions to the imagi- 
nary part of the interaction which was omitted in [2~Tll2"2"] 
and which we keep here to preserve unitarity. 

For the channel with the quantum numbers 1 ++ , we 
are interested only in the S wave in the final state. We 
therefore use the projectors 



1 



1 SS — ^Z ™"'' 1 ' 



mi 
DS 



1 



3n n n n >) (28) 



to find the a^r and af k matrix elements of the ampli- 
tude, where the superscripts S and D denote the S- and 
D-wave components: 

a? k n '(p,pi,E) = a? k s (p,p\E)TS%' + ag s (p,p> , E)T%£ . 

Finally, we calculate the differential production rate 
in the D°D°ir° channel. Further details as well as an 
explicit formula for the production rate can be found 
in [37] ■ A near-threshold singularity (a bound state in 
neglect of imaginary part of the potential), if exists, re- 
veals itself as a below-threshold peak in the production 
rate. 



4. RESULTS AND DISCUSSION 

The aim of the present research is to check whether 
the DD* system interacting through the OPE possesses 



a bound state near threshold. As the first step, we fix 
phenomenologically adequate values Vq = 1.3 MeV and 
A = 1 GeV and find that no bound state exists in the 
system, that we interpret as a demonstration that OPE 
is not sufficiently binding to be responsible for the for- 
mation of the X. To have a better insight, we keep Vo 
fixed and increase the cut-off until a bound state appears 
at the threshold, that is with zero binding energy, that 
gives us the minimal value A m j n = 3800 MeV compatible 
with the existence of a bound state in the DD* system. 
We repeat the same calculation in the static limit, as 
explained before, and in neglect of the OPE imaginary 
part, to arrive at the minimal cut-off A^ = 2750 MeV. 
The difference in the A m ; n between the full treatment 
and the static approach must be ascribed to the effect 
of dynamical pions, inclusion of the imaginary part of 
the OPE potential as well as to the momentum depen- 
dence of the coefficients a. Notice that in the calcula- 
tion performed only the first ordering potential V\ was 
kept. However, large values of the cut-off parameter 
A found indicate that large momenta (with |p| 3> m^) 
may float in the .D-meson loops, and therefore the sec- 
ond ordering potential V2 may give a sizable contribu- 
tion. We therefore supply the OPE potential (fTT|) with 
an extra contribution, coming from the second ordering, 
for which we stick to the static form. Technically this 
amounts to adding to V™ {p,p') from (fTT|) a second 
term of the same form, however with D* and D masses 
interchanged and with the static limit applied, as ex- 
plained above. As a result we find A m i n = 2100 MeV 
and A^ = 1650 MeV, respectively. For convenience 
we collect our results in Table [1] 

One can draw two conclusions from the results ob- 
tained. On one hand, even inclusion of the second or- 
dering term and resorting to the static limit does not 
bring us to the phenomenologically adequate values of 
the cut-off parameter A, the minimal values required for 
the latter being essentially larger than 1 GeV. Besides 
that we observe that the contribution of the second or- 
dering appears to be quite significant both in the static 
limit and in the full treatment. In this situation one 
should, in principle, go beyond the static approxima- 
tion for the second ordering contribution too that, inter 
alia, opens a Pandora box of intermediate states via the 
transition chain 



D*D*n O DD*nn O DDnnn . 



(29) 



undermining in such a way the naive few-body treat- 
ment of the problem. 

Finally, we repeat the same calculations for Vo = 
2.6 MeV which, as explained above, is the value of the 
coupling which was effectively used in 22 . For self- 
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ystat 


^rstat _|_ ystat 


V x 


Vi + v 2 stat 


V =1.3 MeV 


2750 


1650 


3800 


2100 


V =2.6 MeV 


1350 


950 


2400 


1050 



Table 1. The minimal cut-off parameter A m i n (in MeV) consistent with a bound state in the DD* system, that is the 
cut-off that, for a given coupling parameter Vb, provides a bound state residing exactly at the D°D*° threshold. 



consistency, in the full calculation, the values of To and 
r c used in the two-body propagators (|27p were also in- 
creased by a factor of 2. We reproduce the result of 
|22] if, in addition, the static approximation is used, 
imaginary parts of the OPE are neglected, and both 
orderings are taken into account. Namely, a bound 
state at the threshold appears for relatively small val- 
ues of the cut-off, A^f* ~ 1 GeV (see Tabled]). How- 
ever, if the full form for the first ordering potential is 
used, the value of A min = 2400 MeV is required for 
the case of the first ordering potential alone, while in- 
clusion of the second ordering in the static form lowers 
the value of the cut-off to phenomnologically accepted 
one, A m i n = 1050 MeV. This is to be confronted with 
A m in = 2100 MeV obtained for the full treatment using 
the value Vq = 1.3 MeV compatible with the data on 
D* decays. We are led to conclude therefore that one- 
pion exchange is not sufficiently binding to produce the 
AT (3872) state and that other short-range dynamics are 
responsible for its formation. 
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